en 

o 

(N 



X5 



a 



oo 
en 



Liquid-gas phase transition in hot asymmetric 
nuclear matter with density-dependent 
relativistic mean-field models 

Guang-Hua Zhang^ and Wei-Zhou Jiang^'^ 



^ Department of Physics, Southeast University, Nanjing 211189, China 
r~ • ^ National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 



Abstract 



. _ The liquid-gas phase transition in hot asymmetric nuclear matter is studied within 

^ ■ density-dependent relativistic mean-field models where the density dependence is 



introduced according to the Brown-Rho scaling and constrained by available data 



^f) ' at low densities and empirical properties of nuclear matter. The critical temper- 

> ■ 

\^ . ature of the liquid-gas phase transition is obtained to be 15.7 MeV in symmetric 



nuclear matter falling on the lower edge of the small experimental error bars. In hot 
asymmetric matter, the boundary of the phase-coexistence region is found to be 



^— N I sensitive to the density dependence of the symmetry energy. The critical pressure 

^^ and the area of phase-coexistence region increases clearly with the softening of the 

symmetry energy. The critical temperature of hot asymmetric matter separating 
►v> ' the single-phase region from the two-phase region is analyzed to have a moderate 

S I sensitivity to the symmetry energy and is higher for the model possessing the softer 

symmetry energy. 
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1 Introduction 



The determination of the properties of hadronic matter at finite temperature and den- 
sity is a fundamental problem in nuclear physics. At low densities, the so-called liquid- 
gas (LG) phase transition in nuclear matter may occur at sufficiently low temperatures 
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due to the van der Waals behavior of the nucleon-nucleon interaction. The study of 
the LG phase transition in intermediate energy heavy-ion colhsions is of considerable 
interest over the last three decades. In the past, properties of the nuclear LG phase 
transition have been explored both experimentally and theoretically in a variety of 
works [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]. The calculated critical temperature of 
symmetric nuclear matter lies in a wide range, e.g., 13-24 MeV, for various phenomeno- 
logical models [8,9,10,11,12,13]. 

To understand better the features of the LG phase transition in hot asymmetric nuclear 
matter, it is imperative to employ the nuclear equation of state (EOS) of asymmetric 
matter. Recent progress in experiments with radioactive beams provides us a great 
opportunity to constrain the nuclear EOS of asymmetric matter. However, the nuclear 
EOS of asymmetric matter, especially the density dependence of the symmetry energy, 
is still rather poorly known [17,18,19]. The study of the LG phase transition may open 
a possible window to constrain the nuclear EOS of asymmetric matter. In the past, 
the sensitivity of the LG phase transition to the density dependence of the symmetry 
energy was explored with the non-relativistic models in Ref. [15]. Similar work with the 
nonlinear relativistic mean-field (RMF) models can be found in Ref. [16]. In this work, 
we study the thermodynamic properties of the LG phase transition in hot asymmetric 
nuclear matter with density- dependent RMF models. 

The density-dependent RMF models adopted in present work [20,21] feature the chi- 
ral limits at high densities with the in-medium hadron properties according to the 
Brown- Rho (BR) scaling [22,23]. These models achieved satisfactory success in de- 
scribing the nuclear EOS of asymmetric matter, the large mass neutron stars, and 
ground-state properties of finite nuclei [20,21]. It was known that the associated pa- 
rameters that describe the in-medium hadron properties are in agreement with those 
from microscopic calculations [24] or those extracted from recent experimental data 
at low densities [25,26,27]. In addition, the BR scaled NN interaction was succeeded 
in yielding the observed ^^C beta decay suppression [28]. Therefore, it is appealing to 
adopt the models with the BR scahng to study the properties of the LG phase transi- 
tion that occurs in the low-density region. The emphasis is put on the influence of the 
asymmetric nuclear EOS, especially the density dependence of the nuclear symmetry 
energy, on the properties of the LG phase transition. 



2 Formalism 



The Lagrangian density of the density-dependent RMF models is written as 



£ = ^ ^t^.d^ -M*+ gla - gl^y^ - g;i,T,h^,\ V 

~\b,,B^^ + \mfh,,hl (1) 

where ip,a,u, and Bq are the fields of the nucleon, scalar, vector, and isovector-vector 
mesons, with their in-medium scaled masses M*,?7i*,?7i*, and m* respectively. F^^,, 
Bfj^y are the strength tensors of uj and p mesons, respectively 

F^u = d^Ui, — dyUj^, B^y = d^boi, — d^bo^, (2) 

g*, (yf*, and g* are the coupling constants of the scalar, vector and isovector-vector 
mesons with nucleons, respectively. The meson coupling constants and hadron masses 
with asterisks denote the density dependence, given by the BR scaling [20,21,22]. 

In the mean-field approximation, the energy density and pressure at the finite temper- 
ature are given as 



+ E j7^Jd'kE*[n^{k)+n^{k)], (3) 

T=p,n \^'^) '' 

and 



where C^ = gl/m*^, Cp = g;/m*p, C, = nfjgl, E* = \jP + m^2 ^j^h m^ = M* - g^a 
the effective mass of nucleon, a is the isospin asymmetry parameter being a = {pn — 
Pp)/Pb with Pb = Pn + Pp, and Eq is the rearrangement term originating from the 
density dependence of the parameters. 



^. = -p\cfp^-pWCpp^ (5) 

OpB OpB 

OpB ^ ' OpB 

The distribution functions nr{k) and nr{k) for nucleon and antinucleon are given as 



nr{k) = {exp[{E*{k) - Vr)/kBT] + l}-\ 

nr{k) = {exp[{E*{k) + z/.)/A;bT] + 1}-^ (r = n,p), 

where the effective chemical potentials for nucleons read 

i^r = I^T- 9^(^ ± 9*pbo + So, (6) 

with + in ± for neutrons and — for protons. The chemical potentials or effective 
chemical potentials for nucleons can be determined from the nucleon densities 



d^k[nr{k)-nr{k)]. (7) 



i2nY 



With the density- dependent RMF models, we can now study the LG phase transition. 
In this study, we do not consider the finite size effect and thus assume as usual the 
matter is homogeneous. The two phase coexistence is governed by the Gibb's conditions: 

fir{T,p^)=fIr{T,p^), (8) 

p{T,p^)=p{T,p^), (9) 

where L and G stand for the liquid and gas phase, respectively. The Gibbs conditions (8) 
and (9) require the same pressures and chemical potentials for two phases at different 
densities and isospin asymmetries. The stability conditions are given by 

§^) <0,or(§^) >0^ (11) 

oa ) T,p \ da J T,p 



3 Results and discussion 



We first discuss the liquid-gas phase transition in symmetric nuclear matter with the 
density- dependent RMF models SLC and SLCd that only differ in the density depen- 
dence of the symmetry energy [21]. In Fig. 1, we show the pressure of symmetric mat- 
ter as a function of nucleon density at various temperatures. At the low temperature, 
the pressure first increases and then decreases with the increasing density. The p-pB 
isotherms exhibit the form of the two phase coexistence, with a mechanically unstable 
region for each. In both models that share the same EOS of symmetric matter [21], 
there appears an equal critical point at temperature T = 15.7 MeV, where dp/dpB = 0, 
and d'^p/d^pB = 0. This temperature is called the critical temperature beyond which 



symmetric nuclear matter can only be in a single phase. The experimental value of the 
critical temperature extracted for symmetric matter is 16.6 ± 0.9 MeV [10], while our 
prediction is situated on the lower error bar. 
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Fig. 1. The pressure of symmetric nuclear matter versus the nucleon density pB at different 
temperatures in models SLCd and SLC. 

In symmetric matter, the LG phase transition should be of first order, whereas in 
asymmetric matter it was recognized to be of second order [5,14]. In hot asymmetric 
matter, we can nevertheless obtain a temperature, denoted as Tasy^ at a fixed a simi- 
larly according to the relation dp/ dps = d'^p/d'^pB = 0. We show in Fig. 2 the Tasy as 
a function of a in SLC and SLCd. It is found that this temperature is clearly lower for 
the softer symmetry energy that gives higher values in the low-density region where 
the LG phase transition occurs. This is understandable because at given temperatures 
it is easier to access to a larger isospin asymmetry in asymmetric matter with a smaller 
symmetry energy of the model. However, the temperature determined from a point of 
inflection is not the critical temperature for asymmetric matter. In deed, with the con- 
served isospin asymmetry of asymmetric matter the pressure and chemical potentials 
in the LG coexisting phase change throughout the transition [5], in contrast to the 
case of symmetric matter. The calculation of the critical temperature T^ in asymmetric 
matter is thus different, and we will calculate the T^ later on. 

Now, we turn to the discussion on the LG phase transition in asymmetric nuclear 
matter according to the Gibbs conditions (10) and(ll). For a fixed pressure, the two 
solutions of the LG phases to the Gibbs conditions form the edges of a rectangle in 
the proton and neutron chemical potential isobars as a function of isospin asymmetry 
a and can be found by means of the geometrical construction method [5,9]. We thus 
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Fig. 2. (Color online) The temperature T^sy versus the asymmetry parameter a in SLCd and 
SLC. 
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Fig. 3. (Color online) Chemical potential isobars as a function of a at temperature T = 
10 MeV. The curves a, b, c, d, and e correspond to the pressures p = 0.235, 0.175, 0.12, 
0.098, and 0.08 MeVfm~^ in SLCd and SLC, respectively. The rectangle, an example for 
the geometrical construction for p = 0.12MeVfm~^, is used to determine the asymmetry 
parameters and chemical potentials in the LG coexistence phase. 

need firstly to obtain the proton and neutron chemical potential isobars. The chemical 
potentials of the proton and neutron, together with the isospin asymmetry parameter, 
are numerically determined at given pressures and densities. Fig. 3 shows the /i„ and 
fip isobars as a function of a at fixed temperature T = 10 MeV for various pressures. 



The solid and dashed curves are for protons and neutrons, respectively. We see that the 
curves for lower pressures are more complicated than those for higher pressures. Using 
the geometrical construction in the chemical potential isobars, we can then obtain 
two solutions at different isospin asymmetry parameters a for two phases under the 
chemical and mechanical equilibriums. As an example, we plot the rectangle for the 
case of T = lOMeV^ and p = 0.12MeVfm~^ in Fig. 3. The solution with the larger 
a defines the gas phase with the lower density, while the solution with the smaller a 
defines the liquid phase with the higher density. 
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Fig. 4. (Color online) The section of the binodal surface at T = 10 MeV in SLCd and SLC. 
The CP, MA and EC stand for the critical point, the maximal isospin asymmetry and the 
point of equal concentration, respectively. 
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Fig. 5. (Color online) The critical temperature Tc versus the asymmetry parameter a in SLCd 
and SLC. 



The pair of the solutions of the Gibbs conditions found by the geometrical method 
forms the phase-separation boundary at a given pressure in hot asymmetric nuclear 
matter, while all pairs build up the binodal surface. In Fig. 4, we show the section of 



the binodal surface at T=10 MeV. The binodal surface is divided into two branches by 
the critical point (CP) (oc, Pc) and the point of the equal concentration (EC) (a = 0). 
One branch is the high-density (liquid) phase, and the other branch is the low-density 
(gas) phase. We see in Fig. 4 that a much higher CP and a much larger boundary of 
the phase coexistence appear in SLCd, compared with those in SLC. This indicates 
that the section of the binodal surface is very sensitive to the density dependence of 
the symmetry energy, because the unique difference between the SLCd and SLC is that 
the former possesses the much softer symmetry energy than that in the latter. Similar 
sensitivity was found in Refs. [15,16]. Also, we see that the softness of the symmetry 
energy in SLCd makes the maximal asymmetry (MA) a little larger, as compared with 
the MA in SLC. 

As shown in the binodal curve in Fig. 4, the system can only be in a single phase 
for a > Oc at a given temperature. Correspondingly, at a fixed a {a = ac) the given 
temperature plays a role of the critical temperature T^ beyond which the system totally 
cannot fall into the liquid phase at all pressures. In this way, one can obtain the critical 
temperatures at various a. In Fig. 5, we depict Tc as a function of a in SLC and SLCd. 
It is seen that the critical temperature separating a single phase from the coexisting 
liquid-gas phase decreases rapidly with increasing the isospin asymmetry parameter a 
> 0.6. We find from Fig. 5 that T^ is higher for the model with the softer symmetry 
energy (SLCd) that gives higher value of the symmetry energy in the low-density 
region where the LG phase transition occurs. Our finding is in agreement with that in 
Ref. [16]. Interestingly, we find that the sensitivity of Tc to differences in the symmetry 
energy is exactly in contrast to that of Tasy as shown in Fig. 2. As seen in Fig. 5, 
Tc is not very typically sensitive to differences in the symmetry energy. In deed, the 
symmetry energy relies also on the temperature. In order to exhibit the dependence of 
the critical temperature on the symmetry energy, we draw in Fig. 6 the ratios of the 
critical temperature Tc to the symmetry energy, the slope and the curvature modulus 
as a function of the isospin asymmetry. The symmetry energy is here expanded at 
density ph = 0.08/m~^ as 

EsyUPB, T) = EsyUPh, T) + ^Xh + ^Xl + • • • , (12) 

where L and Ksym are the slope and curvature of the symmetry energy at ph = 
O.OSfm~^, and Xh = {pB — ph)/ Ph- The reason for choosing ph = O.OSfm~^ is due 
to the facts that the critical density at the critical point is below but close to 0.08fm~^ 
and that this density may be regarded as an average density over the surface and volume 
of finite nuclei [29]. We see from Fig. 6 that the ratios of the Tc to the symmetry energy 
and its slope are just moderately different for the SLCd and SLC, while the distinction 
is surprisingly large for the curvature. This demonstrates that at the given symmetry 
energy and slope the Tc is moderately different for the SLCd and SLC, consistent with 
the result shown in Fig. 5. Though in the right panel of Fig. 6 the large departure can 
be observed for the Tc at given curvatures, the final role of the curvature is restrained 
by the factor Xh which is quite small in the present case. Moreover, we find that the 



slope and curvature terms in Eq.(12) play an opposite role to that of the zero-order 
term in affecting the critical temperature. As a result, the critical temperature in hot 
asymmetric matter is just moderately model dependent, as shown in Fig. 5. 
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Fig. 6. (Color online) Ratios of the critical temperature Tc to the symmetry energy (left), the 
slope L (middle) and the curvature modulus \k\ (right) as a function of the isospin asymmetry. 
The slope and curvature of the symmetry energy are also calculated at density py^. 

At last, it is interesting to point out that the limiting pressure in the binodal surface 
does not appear in the present models. As pointed out in Refs. [9,13], the limiting 
pressure that is a cutoff of the pressure beyond which the pair of the solution does 
not exist in the geometrical construction is a result of the density dependence of the p 
meson-nucleon vertex. However, this seems to be rather model-dependent. In deed, the 
difference in values of the limiting pressure is large in models applied in Refs. [9,13]. 
In some nonlinear RMF models [16] where the density dependence of the p meson- 
nucleon interaction is actually induced by p-meson effective mass, the limiting pressure 
also does not appear. 



4 SUMMARY 



In summary, we have investigated the effects of the nuclear EOS with the density- 
dependent interactions on the liquid-gas phase transition in symmetric and asymmetric 
nuclear matter with the well-constrained RMF models SLCd and SLC that feature in- 
medium hadron properties according to the BR scaling and can provide satisfactory 
description for properties of finite nuclei and neutron stars. The critical temperature 
of the LG phase transition in symmetric matter is obtained to be 15.7 MeV for both 
models that have the same symmetric part of the nuclear EOS, in nice agreement 
with experimental data. With an analytic continuation from symmetric to asymmetric 
matter, we have obtained a temperature T^sy in hot asymmetric matter according to 
the monotonicity of the pressure and demonstrated that this temperature is lower for 
the softer symmetry energy. It is found in hot asymmetric matter that the boundary 



of the phase-coexistence region is very sensitive to differences in the symmetry energy 
in SLC and SLCd, while the dependence on the EOS of symmetric matter is much 
weaker. A higher pressure and a larger area of the LG phase coexistence region appear 
with a softer symmetry energy. The critical temperature of hot asymmetric matter that 
separates the single-phase region from the two-phase region is found to be higher for 
the model with the softer symmetry energy, while the magnitude of the sensitivity is 
analyzed to be just moderate as compared with that for the boundary of the LG phase 
coexistence. 



Acknowledgement 



We would like to thank Lie- Wen Chen and Chen Wu for useful discussions. The work 
was supported in part by the National Natural Science Foundation of China under 
Grant Nos. 10975033 and 11275048 and the China Jiangsu Provincial Natural Science 
Foundation under Grant No.BK2009261. 



References 

[I] D. Q. Lamb, J. M. Lattimer, C. J. Pethick, D. G. Ravenhall, Phys. Rev. Lett. 41, 1623 
(1978). 

[2] G. F. Bertsch, P. J. Siemens, Phys. Lett. B 126, 9 (1983). 

[3] H. Jaqaman, A. Z. Mekjian, L. Zamick, Phys. Rev. C 29, 2067 (1984). 

[4] H. Q. Song and R. K. Su, Phys. Rev. C 44, 2505 (1991). 

[5] H. Miiller and B. D. Serot, Phys. Rev. C 52, 2072 (1995). 

[6] Y. G. Ma et al. Phys. Lett. B 390, 41 (1997). 

[7] W. L. Qian, R. K. Su, and H. Q. Song, J. Phys. G 28, 379 (2002). 

[8] M. Baldo and L. S. Ferreira, Phys. Rev. C 59, 682(1999). 

[9] W. L. Qian, R. K. Su, and P. Wang, Phys. Lett. B 491, 90 (2000). 

[10] J. B. Natowitz, K. Hagel, Y. Ma, M. Murray, L. Qin, R. Wada, and J. Wang, Phys. Rev. 
Lett. 89,212701(2002). 

[II] Ph. Chomaz, M. Colonna, and J. Randrup, Phys. Rep. 389, 263 (2004). 
[12] A. Rios, NucL Phys. A 845, 58 (2010). 

[13] C. Wu and Z. Ren, Phys. Rev. C 83, 044605 (2011). 



10 



[14 

[15; 

[16 

[17; 

[18 

[19 
[20 

[21 

[22 

[23 

[24 

[25; 

[26 

[27; 

[28 



R. Aguirre, Phys. Rev. C 85, 014318 (2012). 

J. Xu, L. W. Chen, B. A. Li, and H. R. Ma, Phys. Lett. B 650, 348 (2007). 

B. K. Sharma and S. Pal, Phys. Rev. C 81, 064304 (2010). 

B. A. Li, L. W. Chen, and C. M. Ko, Phys. Rep. 464, 113 (2008). 

B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 (2001). 
W. Z. Jiang, B. A. Li, and L. W. Chen, Phys. Lett. B 653, 184 (2007). 
W. Z. Jiang, B. A. Li, and L. W. Chen, Phys. Rev. C76, 054314 (2007). 
G. E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 (1991). 

G. E. Brown, J. W. Holt, C.-H. Lee, and M. Rho, Phys. Rep. 439, 161 (2007). 
X. M. Jin and D. B. Leinweber, Phys. Rev. C 52, 3344 (1995). 

D. Trnka, et. al., Phys. Rev. Lett. 94, 192303 (2005). 
M. Naruki, et. al., Phys. Rev. Lett. 96, 092301 (2006). 

B. Schenke and C. Greiner, Phys. Rev. Lett. 98, 022301 (2007). 

J. W. Holt, G. E. Brown, T. T. S. Kuo, J. D. Holt, and R. Machleidt, Phys. Rev. Lett. 
100, 062501 (2008). 



[29] L. W. Chen, Phys. Rev. C 83, 044308 (2011). 



11 



